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It is shown that the universal critical properties of two recently introduced coupled directed perco- 
lation processes can be described by a single rapidity reversal invariant stochastic reaction-diffusion 
model. It is demonstrated that all renormalizations needed for the calculation of the universal scaling 
behavior near the multicritical point can be gained from the Gribov process (Reggeon field theory). 
Consequently the crossover exponent describing the scaling of the linear coupling parameter is given 
by $ = 1 to all orders of the perturbation expansion. 
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Nonequilibrium processes, their stationary states and 
their phase transitions have been of considerable interest 
in natural science as well as in medicine and sociology 
since many years. Such processes can often be modelled 
by growth and decay of populations with spatially local 
interaction rules. The transition between survival and 
extinction of the population is a nonequilibrium continu- 
ous phase transition phenomenon and is characterized by 
universal scaling laws. For the description of transitions 
in systems that show active and absorbing inactive states, 
percolation models play an outstanding role. Some years 
ago it was conjectured |^,^ that Markovian growth mo- 
dels with one-component order parameters displaying a 
transition into an absorbing state in the absence of con- 
servation laws belong generically to the universality class 
of directed percolation (DP). Besides DP this uni- 
versality class includes e.g. Reggeon field theory (RFT) 

f, the contact process [^,|lO| , certain cellular automata 
and some catalysis models It is relevant to 

a vast range of models in physics, chemistry, biology and 
sociology. 

For the description of universal behavior near a cri- 
tical transition it is often useful to model a universality 
class by mesoscopic stochastic processes involving the or- 
der parameter and other relevant fields. In case of the 
DP class a representation by the Langevin equation for 
the density variable of the so-called Gribov process (the 
stochastic version of the Schlogl model ^) IS appro- 
priate. The name Gribov process was coined by Grass- 
berger who showed that RFT is rather a Markov process 
in disguise than a quantum theory |^ . On the level of a 
formulation of stochastic processes by means of path in- 
tegrals, there is superficially no difference between RFT 
and the Gribov process. However RFT uses creation and 
annihilation operators for particles as the principal fields 
in contrast to the densities and their conjugate response 
fields of the Gribov process. Thus, the latter one is more 
adapted to a mesoscopic description. 

We have been studying the Gribov process for some 
time by renormalized field theoretic methods [|l[. By ge- 
neralization to the multicomponent case we have shown 
recently that the Gribov process with many interacting 



species (directed percolation with many colors) belongs 
to the DP class as well In the multicomponent 

model the evolution of the species is described by a set 
of Langevin equations for density variables {ni{r,t)}. 
These stochastic reaction-diffusion equations are such 
that the condition of the existence of an absorbing va- 
cuum state individually for each species is not destroyed. 
Therefore they are at least bilinear coupled in the den- 
sities. We have shown that multicolored directed perco- 
lation exhibits a new interesting phenomenon: Even if 
the dynamics of the different colors is modelled symmet- 
rically, there exists an instability which leads to a diffe- 
rentiation of the species in their active state under coarse 
graining. The universal properties of this color symme- 
try breaking are described by a multicomponent Gribov 
process with in pairs unidirectional coupled species. 

Very recently further interest on coupled directed per- 
colation processes emerged Tauber et al. have 
introduced a model that couples two species unidirec- 
tional (without feedback) by a linear term in the equa- 
tions of motion. Such a term induces always a nonzero 
mean-density of the second species in a state where the 
first species (the parent) is active. Thus only the inactive 
phase of the first species can be truly absorbing. 

The features of the bilinear and the linear unidirec- 
tional coupled models for two species are combined in 
the Langevin equations 



A"'^ = V'ni - [ti + ■jni]ni + Ci , 



A = V n2- [t2 + '-n2 + fnij ri2 + cmi + C2 , 

(C. (r, t) Q (r', <')> = ^-'g'n^ (r, t) S,,6 (r - r') 5 {t - t') 



(1) 



Here, for simplicity, equal kinetic coefficients A, in- 
traspecies couplings g, and noise parameters g' for each 
species are assumed. As remarked above, the linear cou- 
pling oc a destroys the absorption condition for species 2 
in a state where species 1 is active. Therefore Tauber et 
al. argue that one should introduce further couplings in 
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Eqs. (|l|) because they are now generated through coarse 
graining. This leads to the completed stochastic reaction- 
diffusion equations 

A-'^ = V^m - (n + + Ci , 

^ ^=V-2-(r2 + -n2jn2 

+ - y"-i - /2?^2)?^l + C2 , (2) 

with Gaussian correlations of the Langcvin forces with 
mean zero given by 

(Ci (r, t) Ci (r', t')) - A-15'ni (r, t) S (r - r') S (t - t') , 
(Ci (r, t) C2 ir',t')) = \-'j^n, (r, t) 5 {v ~ r') 5 (t ~ t') , 

(C2 (r, t) C2 (r', i')> = (g'n^ (r, t) + f[n, (r, i)) 

x,5(r-r')(5(t-i') . (3) 

Some remarks are in order here. We consider the cou- 
pled DP processes always in the form of Gribov pro- 
cesses, i.e. in terms of stochastic reaction-diffusion equa- 
tions for density variables. These are in general the 
appropriate mesoscopic variables for the description of 
the percolation-type transition between an absorbing in- 
active phase and an active one as in our case. This 
description is legitimate since creation and annihilation 
processes allowing for single-particle annihilation lead to 
strong spatial clustering. The field theory of Tauber et 
al., which as RFT is based on creation and annihilation 
operators for particles on a lattice, is, up to irrelevancies, 
fully equivalent to our mesoscopic approach. 

In the following we focus on the scaling behavior of 
correlation and response functions in the vicinity of the 
multicritical point where the strongly relevant parame- 
ters Ti and a are small. Critical properties are extracted 
from the mesoscopic model by using renormalized field 
theory and the renormalization group equation in con- 
junction with an er-expansion about the upper critical di- 
mension p^|2^. It is convenient to recast the stochastic 
equations of motion as a dynamic functional [2^-|2^] . The 
Langevin equations (||,^ lead to the following dynamic 
functional in renormalized form 



Jc = j dtd'^xX 



-.§2 

~S2 

/: 



ZX-^dt + ZrT2 - ZxV^ + ^Z„ (S2 - S2) 



2AZ\-^dt - Za(T - Yin - Y2T2 - ZbbS/^ 



1 /' 

^-ZiSi + /2^2S2 - 12^281 - Y^1^2 



Si 



Si 
S2 

(4) 



The Si (r, t) ~ rii (r, t) are the suitable rescaled den- 
sity fields such that g' = g. The Si (r, t) denote the re- 
sponse fields (corresponding to the conjugated auxiliary 



variables of the operator formulation of statistical dy- 
namics by Kawasaki |^ and Martin, Siggia, and Rose 
p?! ) . Correlation and response functions can now be ex- 
pressed as functional averages of monomials of the Si and 
Si with weight exp {— J7c}. The responses are defined with 
respect to additional local particle sources hi (r, i) > 
in the equations of motion (|^). These sources lead to 

the additional terms — {hiSi -\- /i2S2^ in the integrand 

of the dynamic functional Jc As usual, the vari- 

ous Z-factors and counterterms A, Yi are introduced to 
absorb the infinities that are generated in perturbation 
theory by taking the continuum limit. Using dimensional 
regularization and minimal renormalization the different 
counterterms arise as pole expansions in e. Note that for 
a complete renormalization a further gradient term with 
a new dimensionless constant b has to be introduced as 
well as further additive renormalizations that are not con- 
sidered in p^ . 

The scaling by a suitable mesoscopic length and time 
scale, and (A/i^) respectively, leads to Sa ^ Sa ^ 
fj^'i/'^ ^ g ^ f^^ f'^ ^ fj^^f^ where e — A — d. Hence dc = 4 
is the upper critical dimension. A further inspection of 
Eq. (|^) and a glance at the perturbation expansion shows 
that the functional J'c includes all possible relevant cou- 
plings with respect to the Gaussian fixed point as the 
starting point of the perturbation expansion and is thus 
renormalizable. This expansion is dimensionally regu- 
larized and minimal renormalized, whereupon the bare 
(unrenormalized) quantities are related to the renormal- 
ized ones by the scheme 

Si = (Si + AS2) , S2 = ^1/2 (S2 + ^Si) , 

A = Z-^ZxX , 

z-^/^z^^z, 

a = Z-^ (z^a + (Yi + AZ^)ti + {Y2 + AZ^)r2 



Ti — Z^ ^ZrTi , 

5 = Z --Z^'Z^g, Z-^Zbh~2A. 



A = Z^^Z-^'^[Zih - ^ (1 - ^) Zgg ~ 2AZ2f2) , 
f[ = Z-'Z-'/'(z[f[ -A{l-A) Z^g - 2AZ'2fC) , 



f2 = Z-^Z-^'''(z2f2-AZgg 
f2 = Z^'Z-'r2(z'j^_AZn 



(5) 



Dimensionless coupling constants are defined as 



' = G^^A'V-^/^ , (6) 



where = T {1 + e/2) / {inf^'^ . As argued in [||, 
the intraspecies renormalization factors Z, Z\, Z-r, Zg are 
functions of u only and therefore given by the well known 
renormalizations of the one-species model. 

Before we consider the calculation of the various renor- 
malizations to one-loop order we introduce a linear, ho- 
mogeneous transformation between the fields. Let us call 
it the a-transformation: 
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Si 
S2 



Si 
32 



aS2 , 

asi , 



Sl 

32 



Sl 



Note that a corresponding /3-transformation 



Sl 
32 



Sl 
32 



-P~32 

-/3si 



Sl 

h 



Sl 
32 



(7) 



(8) 



is exploited to eliminate S2dtSi from the unrenormalized 
form of the functional J7c- Likewise the scaling transfor- 
mation Si — > 7Si , Si — > 7~^Si is employed to specify g 
as the coupling constant for the nonlinear term and the 
noise term of the intra-species parts of the dynamic func- 
tional (|^). The a-transformation of J7c does not change 
the form of this functional. However it leads to new cou- 
pling constants that we denote with an overbar: 



a = (7 + a{Ti - T2) , 

h= fi+ 2a/2 +a{a-l)g , I2 = h + ag 
f[-2af^ + a{a+l)g , f^^f^-ag 



(9) 



The transformed couplings give rise to transformed renor- 
malization factors which we denote by overbars also: 

A = A , Zb — Zb , Zcf = Za , 
Yi=Yi+a {Zr -Z^) , Y2 = Y2-a {Z^ - Z„) , 
■^1/1 = Zifi + 2aZ2f2 +a{a-l) Z^g , 
Z'Ji = Z'J[ - 2aZy^ +a{a+l) Z^g , 
^2/2 = ^2/2 + aZ^g , Z^/^ = Z'J^ ~ aZ^g . (10) 

Via the transformation formulas the possible fixed 
point values of the coupling constants are transformed as 
well. Thus, the form invariance of the dynamic functional 
under the one-parametric a-transformation explains the 
appearance of the fixed lines found in . Therefore it is 
appropriate to define a-invariant dimensionless coupling 
constants as 



Wi ^ u ' (wi + V2) - V2 , 



(11) 



11)2 — U^^^ {V2 + v'2) ■ 

We define rapidity reversal as 

Sl it) ^ -32 i-t) , Sl (t) ^ i-t) . (12) 

It turns out that under this inversion all unprimed quan- 
tities in (@) change to the primed ones and vice versa, 
and the "temperature" parameters ti,T2 are exchanged. 
Furthermore the renormalization functions are related by 



z,(M,M) ^zliW},{v}) , 
ri(W,M) ^Y2{{v'},{v}) , 



(13) 



whereas Z^,A,Zh are invariant under the exchange 
{v} ^ {v'}. Note that Z, Z\, Zt, Zg do not dependent 
from the sets {v} , {«'} at all. Especially if {v} ~ {v'} 



(a property that is not broken by renormalization), the 
functional J7c is rapidity reversal invariant. Then we have 
= Z>ndYi=y2. 

The somewhat lengthy but simple calculation of all the 
one-loop renormalizations leads to the Wilson functions 
of the renormalization group equation. In part they can 
be gained from the results of |0 in the case 6 = 0. If 
we set u to its fixed point value u* = 2£/3 + O (e^) and 
define as usual f3p = 9p/91n/x|o, where p is any of the 
coupling constants and |o denotes a differentiation hol- 
ding the unrenormalized parameters constant, we obtain 
for the /3-functions of the invariant couplings 

8f3wi = 4 {2wi + W2 ~ u*) W2 + Au*w'i 

+l&u*a - u*b (8wi + IOW2) + 9 {u*bf , 

8/3^j = 4 {2w[ +W2- u*) W2 + 4m*u;i 

+16u*a - u*b {8w[ + IOW2) + 9 {u*bf , 

4/3^2 =8{w2-u*)w2+ 4m* (wi + w[) 
+8u* a ~ 6u* bw2 + i {u* b)'^ , 

8u*Pb = U* (Wi + w[) + {11)2 - U*)W2 

+l&u*a ~ u*b {w2 + u) + {u*bf . (14) 

The function a results from the additive renormalization 
A and mixes the renormalized fields in a correlation or re- 
sponse function under application of the renormalization 
group as 



Anp {si, Sl, S2, S2} = {-as2, 0, 0, -asi} . 



(15) 



Here I?in ^ is the usual linear differential renormalization 
group operator given by 



2?ln^ = Hdfj. + HirTdr + (KctC + KlTi + H2T2) da 



(16) 



with C = 1 - Ix, = 1\ - It, i^a ^ 1\- la and 
Ki = —yi — a. The yi results from the additive renor- 
malizations Yi and 7.. = 91nZ../91n^|o. The function a 
is found to be 

a = —{2 [wi + w[) + 2 {w2 — u)w2/u 

-Abw2 + 3ub^)/l6 ■ (17) 



To one-loop order we have the well known RFT results 
7 = — u/4, "fx — ~u/8, ~ —u/2. The new renormali- 
zations yield 



la = 2 (^"^ ~ "'2) , 

\:{u^/'^Vi+V2v'2 



yi = M 

-6bu^/^V2 - 2bu^/'^v'2 + Sb'^u , 

1/2 — 4:(^U^^^v[ + V2V'2^ 

-Qbu^'^v'2 - 2bu^/^V2 + 36\ . 



(18) 
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In order to determine the fixed-point solutions of Eqs. 
(|l^), P* = 0, we impose the condition a* = b* = 0. 
This yields wl — w'{ — 0, with W2 — (instable) or 
1112 — u* (stable). It can easily be checked that these 
solutions are consistent with the full set of Eqs. ( [1^ ) and 

e). These are the solutions found by Tauber et al. 
. Note that on the fixed lines generated from the 
stable fixed point by the a -transformation two special 
points are found: A rapidity reversal invariant one with 
V2 = V2 = —2vl — —2v'{ — \/vF /2 and a minimally cou- 
pled one with vl = v'{ = V2 = 0, ^2 = The former 
is related to the latter through an a -transformation (|^) 
with a = —1/2. 

Now one has to prove stability of the fixed points 
of the full equations (|l^ ) without using the constraints 



0. A linearization about wi 



6 = and 



either W2 = or W2 ^ u* shows that the flow of b is insta- 
ble for 11)2 — 0, whereas it shows full stability of the fixed 
point for W2 — u* . This vindicates the neglect of a, b and 
the corresponding counterterms A and Zhb in ||l^ at the 
stable fixed point. However without further knowledge 
this statement is only correct in the one-loop calculation 
and could be destroyed in higher loop orders. We will 
show that the stable fixed point is given by w* = w'* = 0, 
W2 — u* to all orders of the loop expansion. As a conse- 
quence the fixed point values a*,b*,yl,y2 are zero, and 
7* = 7^ . This leads to a crossover exponent ^ — 1 where 
$ determines the scaling of cr/|ri|*. 

In the following we will demonstrate that a restricted 
model with only one coupling constant comprises all re- 
levant universal features and especially the stable fixed 
point of all the coupled directed percolation processes. 
We revisit the Langevin equations (^) without the addi- 
tional couplings introduced by Tauber et al. (l^. After 
the suitable rescaling of the densities Ui to the stochastic 
variables Si we now get the renormalized dynamic func- 
tional 



dtd'^x A 



X <; si 

-S2 



Z\-^dt + ZrTi - ZaV2 -f (Si - Si) 

Z\-^dt + Z,T2 - ZxV^ + f Z (.S2 - h) 



S2 -ZaG + Zt]s2 Si 



Si 
Si 
(19) 



The last line of this model functional displays the two 
possible types of couplings of the DP processes, namely 
the linear coupling with constant cr and the bilinear cou- 
pling with constant /. It is easily demonstrated that 
the renormalization is complete by noting that the dy- 
namic functional Jm describes the renormalizable model 
of ||l^ in the case cr = 0. For cr 7^ it remains to de- 
termine the purely multiplicative renormalization of the 
relevant operator S2S1. Explicitly our arguments are as 
follows: For cr = 0, perturbation theory generates only 
the ultraviolet-divergent vertex functions Fg^^j. , Fj-^^.^s. 



with 1 = 1,2, and F^^^^j^sg : where only the latter depends 
on /. All the other possible ultraviolet-divergent vertex 
functions are zero and can only be generated by insertions 
of the composed field S2S1. Such insertions make either 
the vertex functions superficially convergent or produces 
the logarithmically divergent F52.si:(S2Si)- The logarith- 
mic primitive divergency of the latter vertex function is 
eliminated by multiplying Z^ defined via the renormal- 
ized composed field by [s2Si]r = Z^S2Si = Z^Z^^ S2S\. 

We know from our considerations above that the sta- 
ble one-loop order fixed point belongs, up to an a- 
transformation, to the minimal model Jrn ( P^ ) with the 
additional constraint f = g. In it was shown that 
this equality holds at the fixed point in two-loop order 
also. Thus one may suspect that Zf — Zg to all orders 
of the loop expansion provided f = g. We will now show 
that this conjecture is correct. Moreover we show that 
Za- — Zt- 

Of course the minimal model J'm ( |l9| ) is a special case 
of the complete model J7c (^ with /i = /{ = /2 = 
and fi — f- A comparison between both renormalized 
dynamic functionals shows that we can set the counter- 
terms A and Zi,b to zero safely in the minimal case. In 
addition we find Yi = I2 = and Z2 = Zf. After an 
a-transformation (^j^Jl^) we obtain new coupling con- 
stants 



= cr + a (ti - r2) , 
fi = 2af + a{a-l)g , J2 ^ f + ag 
f[=a{a + l)g, = -ag . 

and the additive renormalizations 

Yi ^-¥2 = a {Zr - Z„) , 
A = A = Zbb= Zbb = {) . 



(20) 



(21) 



Now we demand rapidity reversal as an additional sym- 
metry by equating fi = f[. This gives rise to the two 
equations a = —1/2, f = g (there exists a second solu- 
tion with a = 0, f = that corresponds to the bilinear 
uncoupled case and the instable fixed line of [0). In 
this rapidity reversal symmetric case the new coupling 
constants are given by 



/i - /{ = -.9/4 , /2 = g/2 



(22) 



and we have from Eq. dglj) Yi = -Y2 ^{Z^ -_Zr) /2. 
Otherwise we know from Eq. (|l^) that Yi = I2 if ra- 
pidity reversal invariance holds. Therefore both additive 
renormalizations must be zero and from Yi = I2 = we 
obtain the equality 



Za — Zr 



(23) 



Thus, cr and the renormalizc in the same way and 
the crossover exponent $ defined by the scaled variable 
cr/|ri|* is given simply by 



$ = 1 



(24) 
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to all orders of perturbation theory. The other equations 
(0® yield 



Z\ — Z'^ — Zl — Z'2 — Z f — Zg 



(25) 



as one would expect on grounds of consistency. 

The discussion of the fixed-point solutions of the com- 
plete model Jc (H) has shown to first order in an e- 
expansion that the stable fixed-point functional belongs 
to the class of restricted models described by the minimal 
functional Jm up to an a -transformation with the iden- 
tification f — g. Thus the fixed point has a higher sym- 
metry (rapidity reversal) to this order. This higher sym- 
metry is preserved under renormalization. It can easily 
be demonstrated by induction (see appendix B of of 
an analogous case) that rapidity reversal symmetry then 
holds at the fixed point to all orders. Hence all universal 
properties of the coupled DP processes may be calculated 
from the minimal model J7m ( 1^) with identified coupling 
constants f — g and renormalizations Z^ — Z^-, Zf — Zg. 
By virtue of an a-transformation the functional is 
equivalent by to the rapidity reversal invariant functional 



J = J dtd'^x A 



X <; si 



ZX-'dt + ZrTi - ZaV^ + ^Zg (.Si - Si)' 
ZX-'dt + ZrT2 - ZxV^ + \Zg {S2 - I2) 



Si 
Si 



-S2 ZrO -I- Zg- (si - 4S2 -I- 4si - ~Si) Si | (26) 

where u = a ~ (t\ ~ ti) jl. Note again that this model 
solely requires the renormalization constants of the one- 
species Gribov process that are known to second order in 
u. We cite the results of the two-loop calculation 



u 7/2/13 31 35 4\ ^, 

Za = 1 + — + -h — n- +0 {u^\ , 

8e 128£ \e 4 2 3/ ^ ^ ' 

2..2.^.1 + 2^4(I-1)+0M. ,27) 

Computations based on the minimal dynamic func- 
tional Jrn (^) are much easier to perform as calcula- 
tions using the complete model (Q). Thus it may be 
possible to find the equations of state for M2 — 5 (£2) to 
second order and check the assumptions made in [|_7| on 
the reexponentiation of logarithms to yield the new order 
parameter exponent fii of that paper (for a calculation 
of the equation of state for M\ = g{s\) to two-loop order 
see |9|). 

The model J (|2^) describes the coupled DP processes 
near the multicritical point t\ = = (t = 0. What is 



needed for a thorough calculation of fii is a theory that 
comprises the limit cr — > 00. Therefore our considerations 
here do not solve the problem addressed in ||l^, namely 
the determination of the scaling exponent fi^ that arises 
in the scaling (77,2) oc a^^~P'^ (ric — t\)^'^ where species 
1 is in its active phase. Tauber et al. calculate (5i by 
reexponentiation of logarithms. However their approach 
relies on the assumption that simple reexponentiation is 
possible. To derive such scaling properties faithfully one 
indeed has to solve the crossover problem cr ^ 00 which 
induces (possibly!) a new scaling at infinity for the cor- 
relations of species 2. Some features of this crossover 
remind us of the crossover from special to ordinary be- 
havior in the theory of surface transitions js^ , with a cor- 
responding to the surface enhancement c and the species 
1 and 2 corresponding to the bulk and surface respec- 
tively. The crossover problem of interest here is thus as 
yet unsolved. 
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